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Abstract 

We study self-diffusion within a simple hopping model for glassy materials. 
(The model is Bouchaud's model of glasses [J. -P. Bouchaud, J. Physique I 
2 1705 (1992)], as extended to describe rheological properties [P. Sollich, F. 
Lequeux, P. Hebraud and M.E. Gates, Phys. Rev. Lett. 78 2020 (1997)].) 
We investigate the breakdown, near the glass transition, of the (generalized) 
Stokes-Einstein relation between self-diffusion of a tracer particle and the 
(frequency-dependent) viscosity of the system as a whole. This stems from 
the presence of a broad distribution of relaxation times of which different mo- 
ments control diffusion and rheology. We also investigate the effect of flow 
(oscillatory shear) on self-diffusion and show that this causes a finite diffusivity 
in the temperature regime below the glass transition (where this was previ- 
ously zero). At higher temperatures the diffusivity is enhanced by a power law 
frequency dependence that also characterises the rheological response. The 
relevance of these findings to soft glassy materials (foams, emulsions etc.) as 
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well as to conventional glass-forming liquids is discussed. 
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I. INTRODUCTION 



The dynamics of systems close to a glass transition remains a central problem in sta- 
tistical physics Because the glass transition is a many-body phenomenon, models for 
it invariably involve some approximation and it is important to disentangle the physical 
phenomena from the approximation scheme used. In mode-coupling theories, for example, 
the dynamics are dominated by a small number of collective degrees of freedom [Q. This 
is appealing, but ignores activated processes, and obscures the fact that glassy materials 
are dynamically heterogeneous with a variety of local environments . Another approach 
is to use disordered hopping models, which treat instead single-particle degrees of freedom 
coupled to a random environment, and this is the route followed here. Such models, though 
unrealistic in some respects, are likely to give a better account of those properties, such as 
self- diffusion, which are not dominated by relaxation of collective modes. 

An important breakthrough in the approach based on hopping models was that of 
Bouchaud [Q, who showed that an ensemble of non-interacting particles, moving by ther- 
mally activated hopping at temperature T in an uncorrelated fashion through a sequence 
of traps, can show a glass transition. This occurs if and only if the density of states (the 
prior probability distribution of trap depths) p{E) has an exponential tail. Though grossly 
over-simplified, this allows many properties of the model to be calculated exactly. The expo- 
nential tail to the prior distribution p{E), which leads to a power-law spectrum of relaxation 
times and is vital to the appearance of a glass transition, is supported by evidence from the- 
ory on spin glasses |^ and experiments on low-temperature fluids [H. Such a power-law 
relaxation spectrum was the point of departure for a related model of glasses 0J^, and 
similar features appear in models of dispersive transport in disordered semiconductors 0. 

The original model of Bouchaud has no explicit spatial coordinates. Nevertheless, the 
extension to self-diffusion is unambiguous (if each hop corresponds to a spatial displacement 
with a well-defined second moment) and indeed Monthus and Bouchaud gave an expres- 
sion for the distribution of displacements on a hypercubic lattice. For flow properties, more 
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choices are possible in the way spatial coordinates are treated, but a minimal extension of 
the model was offered by Sollich et al. |Tl||T^. This model, which we call the GR (glassy 



rheology) model allows both linear and nonlinear rheological properties to be calculated. 
Details of it are recalled in section |T| below. 

The model introduced in Ref. |TT]] was in fact proposed to describe a class of "soft glassy 
materials" (argued to include foams, dense emulsions, etc.), in which context it was found 
necessary to replace the thermodynamic temperature T in Bouchaud's model by an effective 
(noise) temperature x ^ T. This replacement converts the GR model studied below into 
the soft glassy rheology model of Refs. |[TT|Jl^ , and can be made throughout our calculations. 

The glass transition in the GR model shows interesting features which may need fur- 
ther explanation for readers whose background lies in conventional glasses. Specifically, in 
Bouchaud's model the glass transition Tg is identified as the temperature below which the 
system shows "weak ergodicity breaking": its Boltzmann distribution is not normalizable. 
This means the system will evolve into deeper and deeper traps as time goes by, and will 
never attain a steady state — although at no finite time are there infinite barriers parti- 
tioning phase space. When flow-related degrees of freedom are included in the model (see 
section |ITT| below), one finds, as expected, that the system has a finite elastic modulus (at 
zero frequency) for temperatures T < Tg. Less obviously, the viscosity of the material di- 
verges, not at T = Tg but at T = 2Tg. Between these two temperatures, the static modulus 
is zero but the viscosity infinite; the material is what is known in rheological language as 
a "power-law fluid". For applications of the model to soft glassy materials this is a very 
attractive feature, since such viscoelastic behaviour is frequently observed in these systems 



rT|,|T2|. In a model of conventional glasses, on the other hand, it might be considered un- 
desirable. However, with the introduction of a high energy cutoff in the prior distribution 
p{E), the phenomenology of the GR model can be adapted to model that of conventional 



glasses. We discuss the conceptual features of this modification at the end of section |V11|; for 



the actual calculations in the present paper we restrict ourselves to the simpler case without 
cutoff. 



Below we calculate the statistics (additional to the results of Ref. [|T0|) of self- diffusion 
in the GR model, and explore quantitatively the breakdown, near the glass transition, of 
the (generalized) Stokes-Einstein relation [(G) SER] between self-diffusion of a representative 
particle in the fluid and the (frequency-dependent) viscosity of the system as a whole. (This 
breakdown does not, of course, imply that GSER also fails for macroscopic probe particles 
which see the surrounding material as a continuum.) We also investigate the effect of shear 
flow on self-diffusion and show that this causes a finite diffusivity in the low-temperature 
regime where this was previously zero. 

In the next section we review the status of the GSER in both conventional and soft 



glassy systems. The GR model is described in section |T| and used in section |V|to study and 
discuss the breakdown of the GSER. In section 0, we further elucidate the physics of particle 
transport in glassy systems, using both analysis and simulation of the model. In section ^ 
we calculate the effects of shear on the rate of self-diffusion, finding the diffusion constant 
in the presence of continuous and periodic shear strains of various rates and amplitudes. In 
section |Vli| we conclude with a discussion of our results in the context of both conventional 
super-cooled/glass- forming liquids and soft glassy materials. 

II. STOKES-EINSTEIN RELATION 



The Stokes-Einstein relation (SER) between coefficients of diffusion and viscosity (and its 
generalization to frequency- dependent quantities) is of great importance to our understand- 
ing of transport in fluids. This is especially true since the development of new techniques 
to measure viscoelastic response by light scattering from (or tracking of) probe particles 
r3HT5[. Yet the relation's range of validity has been called into question by experiments on 



super-cooled and glass- forming liquids [p!7|-p2|. 

To find the linear storage and loss moduli, G'{uj) and G"{uj), of a viscoelastic fiuid, 
a rheometer is normally used to apply an oscillatory stress to the fiuid and measure the 
resultant strain (or vice versa) as a function of (angular) frequency u. A technique has 
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recently been developed whereby the complex viscoelastic modulus G* (= C + iG") may be 
foundQ, even when the available sample of fluid is too small to flU a rheometer. The technique 
involves introducing one or more neutrally buoyant probe spheres of radius a into the fluid 
and measuring their Brownian motion. For a Newtonian fluid, which is characterized by 
a single, frequency-independent viscosity 77, that viscosity is given by the familiar Stokes- 
Einstein formula rj = ksT /QitaD . Here ks is Boltzmann's constant and T temperature. As 
the fluid is Newtonian, the probe executes an unbiased random walk, and its mean-square 
displacement as a function of time is = QDt^ which deflnes the diffusion constant 

D appearing in the Stokes-Einstein formula. In this purely viscous fluid, G* is purely 
imaginary, given by G*{uj) = iujrj. As the derivative of the stress relaxation modulus is a 
well-behaved function of time, its Laplace transform G{s) is simply related to its Fourier 
transform G*{uj) = G{iuj). Combining the above expressions, we see that probe spheres in 
a Newtonian fluid respect the relation 

nas {r [s)) 

where (r'^{s)'j is the Laplace transform of the mean-square displacement of a probe. Equation 
D says that the modulus of a viscoelastic fluid is the ratio of a driving force (in this case 
induced by ksT rather than by a rheometer) to a response (mean-square distance moved by 
the sphere). 

Although derived above only for the Newtonian case, Eq. |l] is the generalized Stokes- 
Einstein relation (GSER), and when written in this form it applies under relatively general 
conditions^ to all linearly viscoelastic fluids, not just Newtonian ones |T^. The basic Stokes- 
Einstein relation (SER) is recovered from Eq. |I| in the zero frequency limit, since viscosity 

^G*{u)) is the Fourier transform of the time derivative of the stress relaxation modulus Gr{t), 
which relates stress on the fluid, a{t) = /g dtiGr{t — ti)j{ti), to strain rate 7(t). 

^At sufficiently high frequency, for a probe of finite mass, the modulus has an additional inertial 
contribution which we drop. 
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is defined by 77 = Mm^^o G" (u) / u . In all cases, the fluid around the probe is treated as a 
featureless continuum, and therefore the result is only valid when the sphere radius greatly 
exceeds the microscopic length scale {e.g. particle size) of the fluid. (Furthermore, probe 
particles must be sufficiently dilute to ensure independent Stokesian flow fields around each.) 
If, as is sometimes the case, the GSER is applied to probe particles which are not large com- 
pared to the microscopic scale, then it ceases to be a rigorous result. As an approximation, 
it is then close to the spirit of effective medium theory (wherein a representative particle is 
viewed as embedded in a continuum and its contribution to the continuum properties then 
calculated self-consistently). 

With macroscopic probe particles, the GSER has been fruitfully employed p3|JT^ to mea- 
sure the rheological properties of viscoelastic fluids from light scattering or optical tracking 
observations of the suspended probes. In Ref. |15|, Mason and Weitz applied this method 
successfully to various systems, using colloidal spheres as probes. Remarkably, in one case 
the medium analysed was itself a suspension of the colloidal particles, of volume fraction 
= 0.56, which lies close to the glass transition concentration. The scattering method 
used was diffusion-wave spectroscopy (DWS) |^ whose direct interpretation in terms of 
self-diffusion is somewhat ambiguous; but in any case, both criteria for the rigorous validity 
of the GSER were violated (no separation of length scales, and probe particles not dilute). 
Yet the diffusion data for this system, analysed as though DWS was probing self-diffusion 
and transformed into putative rheological data via the GSER, compare remarkably well with 
frequency-dependent rheological data measured directly. Among various explanations that 
are possible for this, is the idea that the GSER is in fact more widely valid than generally 
assumed ||T5[| . 

On the other hand, in conventional glass-forming fluids, the (zero-frequency) tracer dif- 
fusion and viscosity have been measured as a function of temperature, by experiment [p!^^-|22 



and simulation |^3|]. While the SER held (within a factor of order unity) for probes just a few 
times larger than the fluid particles |2D|, for equal-size probes the translational diffusion was 
enhanced by orders of magnitude with respect to the SER prediction as the glass transition 



was approached. 

In view of the unexpected success of GSER for colloidal glassy fluids, alongside the failure 
of SER for conventional glasses^, a more detailed examination of the connection between 
diffusion and rheology seems worthwhile, and we pursue this in what follows. 

III. THE MODEL 

We use a simple hopping model to emulate properties of a fluid close to its glass transition. 
The model, which will be referred to as the 'glassy rheology' (GR) model, was defined in 



Ref. 1 11], where its response to a macroscopic shear was calculated. In the present study we 



shall investigate the relation between rheology and particle transport in the GR model. 

A. Definition 

The GR model is based on Bouchaud's model for glassy dynamics ||^, with the addition 
of shear strain degrees of freedom which that model lacks. Particles in the Bouchaud model 
are thermally activated at temperature T = (3~^ from traps of energetic depth ii^ in a time r 
which has an exponential distribution with rate Fq exp — {(3E), where Tq is some microscopic 
rate constant. Having escaped, a particle selects its next trap at random {i.e. there are no 
spatial correlations) from the prior distribution of trap depths p{E). We briefly recall the 
resulting behaviour of Bouchaud's model. For high temperatures T, the occupancy of traps 
of depth E evolves towards the Boltzmann distribution P^^^{E) ~ p{E) exp{E/T). As T 
is lowered, this distribution may cease to be normalizable, leading to a glass transition at 
T^^ = — \imE->oo{d/ dE) In p{E). For T < Tg, no equilibrium state exists, and the system 
shows "weak ergodicity breaking" and various ageing phenomena. A finite value of Tg implies 
an exponential tail in the density of states, p ~ exp{—E/Tg). 



^For a related discussion of the failure of Einstein relations in biased diffusion processes, see |24| 
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The GR model ascribes some internal features to the traps. This is done through a 
minimalist description in which the tensorial aspects of the problem are ignored and the 
shear rate and strain are treated as though they were scalar quantities |]Tl],|12|. The trap's 
potential is taken to be quadratic in a local strain (or relative displacement) coordinate I, 
so that the particle's energy is ^kP, with k a constant. This applies up to a local yield 
value /y, whose energy E = \kly^ is the depth of the trap. Upon reaching this threshold, the 
system locally rearranges to a new configuration, thus relaxing the local strain ly. The yield 
energy of the new potential well is drawn from the prior distribution p{E), and its strain / 
set to zero. In the absence of an externally imposed shear rate 7, all I values remain zero, 
and the GR model reduces to the Bouchaud model, with hopping between traps caused by 
thermal activation alone. With shear imposed, changes in all local strains are assumed to 
follow the global strain rate, / = 7. (This assumption is mean-field in character.) With 
the particles thus dragged up their quadratic energy wells by global shear, the barrier to 
thermal activation is reduced, and the local yield rate becomes roexp[— /^(i? — ^kP)]. If 
Pocc{l, E; t) dl dE is the probability of finding an occupied trap of given I and E, the above 
dynamics imply [ p!T] , p!2| 



^Pocc = -T^^cc - V.e-^^'^-^'"^ Pocc + T(t) p{E)S{l) (2) 

The first term on the R.H.S. arises from local straining by the macroscopic deformation. 
The second describes the loss of occupied traps from the distribution by thermally activated 
yielding. The last term corresponds to creation of new, unstrained configurations, whose 
trap depth is drawn from the prior p{E), at a rate equal to the total yield rate in the system, 
r(t)=ro {exp[-f3iE-lkP)]] 

Note that, as originally conceived in the context of soft glassy materials each 



occupied trap represented a small region of the medium, large enough for a local shear 
strain / to be defined, but small enough for this to be approximately uniform within the 
region. However, the model can readily be interpreted as describing the motion of individual 
particles in a conventional super-cooled/glass- forming fluid, by associating the local strain 
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/ with, say, the relative displacement of a particle from the centre of a trap formed by its 
neighbours. 

We have not yet defined the stress in the GR model. Due to stochastic yielding events, the 
local stress, which we take to be kl, is inhomogeneous (it follows the local strain /, rather than 
the macroscopic strain). Since we have stated that the strain rate / is everywhere equal to the 
global strain rate 7, these local stresses must combine according to a straightforward average 



so that a = k {I) = k J I Pocc{l, E; t) dl dE. As discussed in section p^V A| , this corresponds to 
a fully parallel mechanical circuit; this is the only combination consistent with affine shear, 
and may well be approached in many physical systems. 

B. Rheology and Diffusion 



The dynamics described (Eq. ^ are further motivated in Refs. |TT|JT^, where the GR 
model's constitutive relation between stress and strain rate is calculated. We now briefiy 
summarize the linear rheology that the model predicts. We use non-dimensional units for 
time and energy by setting Fq = Tg = 1; we also re-scale our strain variables (/, 7) and stress 
(J so that k = 1. The complex dynamic shear modulus G*{ijj) = G' + iG" describes the 
stress response to small shear strain perturbations around the equilibrium state. As such, 
it is well defined {i.e. time-independent) only above the glass transition, T > 1. Expanding 
Eq. (H) to first order in the amplitude 70 of an oscillatory strain 7(t) = 70 cos tut, we find 

= / ^"^^ ) (3) 

where the average is taken with respect to occupied traps, and t{E) = exp{f3E) is the mean 
(Arrhenius) residence time for a trap of depth E. (We distinguish this from the actual 
residence time r in such a trap, which is a random variable.) Equation |^ corresponds to a 
distribution of Maxwell modes whose spectrum of relaxation times is given by the equilibrium 
distribution P^^^{E) ~ exp{PE)p{E). Given that p{E) has an exponential tail (required for 
Tg to be finite), the distribution of residence times thus exhibits power-law behaviour for 
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large r: \1/occ(t) ~ r ^. This leads to power lawsQfor G* in the low frequency range: 
G' ~ uj^ for 3 < T, ~ uo^-^ for 1 < T < 3 

(4) 

G" ~ for 2 < T, ~ cj^-i for 1 < T < 2 

For T > 3 the system is Maxwell-like at low frequencies, whereas for 2 < T < 3 there is 
an anomalous power law in the elastic modulus. Most interesting is the regime 1 < T < 2, 
where G' and G" have constant ratio; both vary as uj'^~^. Moreover, the frequency exponent 
approaches zero as T ^ 1, resulting in essentially constant values of G" and G' . Note, 
however, that the ratio G" / G' ~ T — 1 becomes small as the glass transition is approached. 
This increasing dominance of the elastic response G' prefigures the onset of a yield stress 



for T < 1 11.12 . However, for T < 1 the linear viscoelastic moduli show slow time evolution 



and ageing effects |2^ in accordance with the weak ergodicity breaking discussed in section 
|. Accordingly, in the following section, where the GSER for this model is discussed, we 
consider only T > Tg = 1. However, that does include the regime Tg < T < 2Tg which 
would lie below the "glass transition temperature" were this to be defined operationally 
by the divergence of the viscosity. As Tg is approached from above, the equilibration time 
following a quick quench diverges. Hence, even for T > Tg, non-equilibrium situations 
are also of interest in glass-forming materials, though we shall consider only equilibrium 
situations in the following section. 

In order to study self-diffusion we must clarify further the spatial interpretation of 
our GR model. We shall associate each relaxation event with a step on an unbiased random 
walk. We draw each step from a GaussianQ of unit variance in each spatial direction (thus 



^Such phenomenology is not confirmed for conventional glass formers. However, with a high- 
energy cutoff in p{E), which may be apparent in conventional glasses, Maxwellian behaviour re- 



appears at low uj. See section VII 



^Given that the walk is Markovian, our choice of geometry of the individual steps will not influence 
the large-scale behaviour. 



11 



defining a length unit for tlie model). This gives a meaning to the square displacement, r^(t). 
After each step, the randomly walking particle is temporarily trapped within a potential well. 
Within each well, we assume rapid 'thermalization' of the particle, which therefore samples 
a Gaussian distribution of displacements at each site on the walkQ. In fact, we assume this 
thermalization to be instantaneous, thus discarding information about high frequency intra- 
well dynamics. Hence, we model the time regime of a-relaxation, but not of /^-relaxation. 



IV. TEST OF THE GSER FOR SMALL PROBES 

We have the complex linear viscoelastic modulus for the GR model, quoted in the section 
above. It is now our task to compare this with the prediction of the GSER, applied to a 
probe particle which is representative of the fluid itself. The result of Eq. 0, whose validity 
is in doubt for a small probe of this kind, will be denoted Ggp- Calculating Ggp from Eq. 
|1] requires a knowledge of (r'^{s)'^ in the equilibrium fluid. The brackets (...) denote an 
average over all random walks of the particles, in terms of choice of path and of trap depths 
drawn from the prior distribution p{E). In the absence of intra- well structure, the treat- 
ment of continuous-time random walks is discussed in Refs. p4| , p6| , ^ , p!0| . The mean-square 
displacement can be found from the correlation function C, defined by C{q,t) = {^e*^'^*^*^ 
This is calculable 0] via C{q,t) = J2N=oQNiQ) Pt{N) where Pt{N) is the probability of 
having performed exactly hops at time t, and Qn{q) is the Fourier transform of the po- 



sitional probability distribution after an A^-step random walk. As described in section |T|, 
the positional distribution is the convolution of a standard random walk and a "Gaussian 
of thermalization" (whose variance is Td/k, where = 1 is the local spring constant, and d 
the dimension of space). So its Fourier transform Qn{<i) is the product of the usual function 



^We neglect modification of the Gaussian by the leaky boundary condition at the edge of the 
well. This is consistent with the unmodified hopping probability used to construct the GR model 
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exp(— and another Gaussian. 

The quantity Pt{N) depends only on the distribution of residence times. Note that the 
time between hops along a random walk is drawn from the distribution \E'hop(T), which is 
not the distribution of residence times of occupied states in the system, \&occ(t). The former 
distribution gives the inter-hop times available to a random walker upon selecting its next 
trap. The latter gives the residence times of all the states in the system which are presently 
occupied. Thus \E'occ(t) is related to ^^hopl^) by a weight factor r, since the hkelihood of 
finding a given trap occupied is proportional to the time for which it is occupied: 

= i^:^. (5) 

(^)hop 

Accordingly, averages of some stochastic quantity Q with respect to the two distributions 
are related by 



(Q)hop = i^^, (QU = ^77r^. (6) 

In terms of a typical random walk, (Q)hop is the average with respect to the steps, while 
{Q)occ i^ average with respect to time. Note that, in calculating the probabihty Pt{N) 
of having performed hops at time t, the time spent in each trap is drawn from \l'hop(''") 
except for the first trap |^^, which is selected from \E'occ(t)- This reflects the fact that 
the walker is not introduced to the system at time zero, but is selected at random from the 
traps already occupied. This careful choice of the first trap is important to the diffusive 



behaviour at early times and, as we shall see, its influence persists for increasingly long 
times as the glass transition is approached. 

From the definition of C{q,t) it follows that, in terms of its temporal Laplace transform 
C{q,s), 



r'{s)) = -VijC{q,s)^^Q (7) 
where Vq denotes a derivative in qr-space. From this prescription, we find 



r'^is' 



s(l-^hop(s)) s 
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Here, "^{s) is the Laplace transform of ^I/(t). The distribution of available escape times 
^hop('i") is given by averaging the exponential distribution for a single trap over the prior 
energy distribution p{E) thus: 

POO 

^hop(r) = / dEp{E) exp -[(3E + re"'^^]. (9) 
Jo 

So its Laplace transform is 

with t{E) = e^^ the mean residence time for trap depth E. The expression for \l'occ('S) is 
the same, with the average over hops replaced by an average over occupied traps. So, from 
Eq. I, 

$„,a.)^^P^f<^^i^ffW (11) 

(^)hop ^0 1 + sr{E) s (r)h„p 
Substituting this expression into Eq. ^ cancels the \Ef-dependence. So we see that the second 

moment of the distribution of displacements, 

is independent of any properties of the distribution of trapping times except for its mean. 

If we use this expression for the self-diffusion of a small probe to calculate the modulus 
in Eq. |l| (and substitute s = iuj to find the complex viscoelastic modulus in the Fourier 
domain), we findQ that the GSER, applied to a small probe, predicts a complex modulus 



With ro = T (r),„p = T/ (r-\ 

This is exactly Maxwellian behaviour at all temperatures for which an equilibrium distribu- 
tion exists {i.e. for all T > Tg = 1). Clearly, this expression is completely at odds with the 



actual rheology of the GR model as described in section III. 



^We have set d = 3 and particle radius a = kB/Sn to simplify the constant of proportionality in 
the GSER (Eq. 1). 
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A. Breakdown of the GSER 



The generalized Stokes-Einstein relation, when applied to a small probe, has failed in 
two distinct ways which we now discuss. The first is a (fairly trivial) discrepancy in the 
temperature dependence, which is inherent in any hopping model. To study the spurious 
T-dependence of the characteristic time in Eq. |T3|, let us take the inverse Laplace transform 
of Eq. 0. Thus we find the nature of the probe's diffusive motion as a function of time: 

{r\t)) = Td + td/{T)^^^. (14) 

The time-independent term Td is the variance of the local 'Gaussian of thermalization' in 
which the probe finds itself at each step of its walk, and respects equipartition of energy 
amongst the elastic degrees of freedom. The second term says that diffusion is exactly linear 
on all time scales^ with a diffusion constant 

The GSER implies that the mean-square speed (also given by equipartition, (||mf^^ = 
^kBTd) sets the rate of translation along the random walk. Indeed this is why the diffusion 
constant is proportional to temperature in dilute gases, where particles are not caged, so 



their directions rapidly become decorrelated. Like any activated hopping process ||29| with 
a temperature-independent microscopic attempt rate {e.g. the Eyring model of dense fluids 
poll), the GR model lacks such a factor T in the diffusion constant (Eq. |T5]). This is the 
source of the rogue T-dependence of G*p. As the GR model is designed to be applied close 
to the glass transition, factors of T are of order unity. This source of deviation from the 
GSER is therefore not of great importance to most experiments, which usually only test it 
up to factors of order unity |^|T3|,|TB| . 



^We note that the diffusive term in Eq. is equivalent to the infinite-dimensional limit of the 
exact result for a random trap model on a hyper-cubic lattice obtained by Schroeder [^] and Kehr 
et al. [|6|,p. 
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However, the GSER is violated by much larger factors than this arising from the second 
source of error in Eq. |T3|. This is the "misplacement" of the thermal averaging brackets, when 
compared with the correct expression in Eq. ^. In a Maxwellian fluid, the distribution of 
lifetimes is narrow, since there is a single characteristic relaxation time. Thermal averaging 
then commutes with the other algebraic operations, so the incorrect placement of the thermal 
brackets becomes irrelevant in that case^. However, the GSER result breaks down for broad 
distributions such as exist both in many real glass formers, and in the GR model. To 
illustrate this, let us compare the diffusion constant (Eq. |15]) with that predicted by the 
Stokes-Einstein relation (zero- frequency limit of the GSER), given the GR model's actual 
viscosity. From Eq. ^ with t] = lim.i^^Q G" (u) / u , the SER gives the diffusion constant 
as T/2t] = T/2 (t)^^^. The first moment of the distribution of times for which traps are 



occupied {t)^^^ diverges in the GR model for temperatures below T = 2, whereas (r ) 



occ ' 



which appears in Eqs. |T5| and remains finite. Thus D is enhanced with respect to the 
SER value as the viscous divergence is approached, as is widely observed in experiment and 
simulation |]T^-^,^]. We note, however, that the GR model is not sufficiently elaborate to 



account for the power law {1/D oc r]^,C, < 1) which has been observed to replace the SER 
in some supercooled fluids . 



The rheology of the GR model is equivalent to a set of over-damped harmonic oscillators 
[e.g. masses on springs in dash-pots) connected in parallel. Each over-damped element has 
the characteristics of a Maxwellian fluid, with a different time constant. Even if only a 
vanishingly small fraction of the population of Maxwell models becomes infinitely viscous 
{i.e. their time constants r diverge), it will be apparent in the response of the system, 
whose stress is the sum of stresses of the population^. Contrast this with the GSER for a 



^In fact, the GSER for small probes is tantamount to mean field theory, in which averaging indeed 
commutes with certain algebraic operations. 
^°The connectivity of the damped oscillators, which is parallel for the GR model (section |l|), 
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small probe, Eq. 0. It claims to derive the moduli from a knowledge of only the second 
moment, (r^), of the distribution of displacements. This is (qualitatively) more like summing 
compliances of the population^ 

This discussion shows that the second moment alone, (r^), does not contain sufficient 
information about the distribution of displacements (or, equivalently, about the full distri- 
butions ^^(t)) to find the rheological behaviour of the GR model. Specifically, this moment 
would be unaffected by a vanishingly small subset of the population being stuck for infinite 
time at r = 0; yet this subset may dominate the rheology, e.g. by causing {t)^^^ to diverge. 

V. HOPPING STATISTICS IN THE GR MODEL 

As we have mentioned, the GR model shows some intriguing dynamical features. The 
viscosity diverges at a temperature (T = 2) well above the glass transition temperature 
[Tg = 1 defined as the temperature below which the system has no equilibrium steady 
state). The diffusion constant D vanishes at Tg, and yet not all particles are static below 
this temperature. In this section we gain a fuller picture of the model's dynamics, by 
investigating its hopping statistics analytically and by simulation^. 

Let us begin by comparing the two distributions \l/hop(T) and \1/occ(t), for the residence 
times of particles in a trap. The first is taken over the prior distribution (of a priori available 



will be different in many real systems, becoming a more elaborate network of series and parallel 
connections. 

^^Of course, the valid application of the GSER (to a large probe) is equivalent to parallel con- 
nection of Maxwell models, since the surface of the large probe feels the simultaneous influence of 
many degrees of freedom, which, within the GR model, have additive stresses. 

^^A similar analysis was carried out by Kehr et al. |26,2^ on a real lattice {i.e. with spatial 
correlations). They modelled diffusion in crystals which exhibit no glass transition, and therefore 
used a distribution of residence times narrower than ours. 
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states); the second over those states which are actually occupied at any given time. As stated 
in section |IV|, these two distributions are so different that selecting just the first step on a 
random walk from the wrong distribution can have drastic effects on {r^{t)) which persist up 
to rather late times. Suppose that by mistake the initial residence times were drawn from 
the distribution \I'hop(T) (which is the correct choice for all subsequent residence times), 
rather than from \1/occ(t) (which is correct for the first one only ||26|j27[| ). The mean square 



displacement ('" (i)) prior '^^^^ ^^^^ hypothetical ensemble has Laplace transform 



7^{s)) + ^ (16) 

Vprior ,(1-Vl>,„p(s)) S ^ ^ 



which, in contrast to Eq. |T2|, does not depend only on the mean of the distribution \l/hop(T)- 



The resultant (?"^(t)) prior ^ non-linear function of time (tending to Eq. [T^ as t — oo p4| ) 
which increases monotonically, even at T = Tg (where the true diffusion constant vanishes). 

The correct value of (r^(t)) given in Eq. ^ is not found as just outlined, but instead 
by weighting each walk in the ensemble with the residence time of the trap it is initially 
in |^,27|. This factor is required because the probability of a particle initially being found 



in a given trap is proportional to the residence time of that trap. (Since we are interested 
in a system which has already equilibrated, the particles have, in principle, had time to 
sample the distribution of traps, and are therefore more likely to be found in long-lived 
states than a priori.) It is somewhat unintuitive (though necessary as a consequence of the 
system's steady state) that this re-weighting of the ensemble should conspire to cancel all 
non-linearity from the time dependence of mean square displacements implicit in Eq. |l^ and 
yield Eq. instead. 

For the exponential prior p{E) = exp{—E), it can be deduced straightforwardly from 
Eqs. |10| and ^ that, in Laplace time. 



^ poo 

^hop(s) =T du 
Jo 



s + e-" 

Tu 



^occ(5) = (T-l) / du^^ 

Jo s + e ^ 

from which it follows, showing temperature dependence explicitly, that 



^occ(r; T) = ^hop(r; T-1). (17) 

Thus, at equilibrium, the residence time distribution of occupied states coincides with the 
prior distribution of such times at a substantially lower temperature. 

A. Numerical results 

To observe the process in action, we simulated the activated hopping of 10^ particles, 
using the exponential prior distribution of trap depths p{E) = exp{—E). For each particle, 
the time interval for successive hops was drawn from \l/hop(T) (Eq- 1^) for all but the first hop 



(as in p6[). The algorithm used to simulate hopping in the system between times and t is 
as follows. Each of the 10^ particles is initially scheduled with a hopping time drawn from 
^occ(t). If this is greater than t, the particle is not visited again, and is recorded as per- 
forming zero hops during the simulation. The remaining particles are successively assigned 
further relaxation times[^ drawn from \l/hop(''") and their numbers of hops incremented, until 
the next hop is scheduled for after t. In this way the distribution pt{N) of the number of 
relaxation events up to time t was measured, as a function of temperature. 

The distribution of hops Pt{N) gives a clearer picture of dynamics in the system than the 
distribution of displacements, which is just its convolution with a random walk for which 
(r^) = Nd. The results appear in Fig. |l| for a range of temperatures and times. We see that, 
at high temperature, all the particles are mobile and concentrated in a peak whose position 
increases linearly with time. As the temperature is lowered towards the glass transition 
(T = 1), a second peak containing a significant fraction of the population appears at = 0. 
Of course, as this peak decays, it feeds the rest of the distribution. Its lifetime increases 
as the transition is approached, yet the rest of the distribution, for which it is the source, 
remains concentrated in a mobile peak. 



^^The stochastic variable t = — ^ ^/-^In,^' is distributed according to ^'hop(7") if ^ and ^' are 
distributed uniformly on (0, 1). 
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B. Analysis of bimodal diffusion 



We consider next the time dependence of PtiO), which is the probabihty that the time of 
the first hop exceeds t: 

Pt{0)= ^occ{T)dT. (18) 

Using Eq. and expressing \l/hop(T) in terms of an incomplete gamma function, \l/hop(T) = 
r-(^+i)7(T + 1, r)/T (from Eq. | with p{E) = e'^), we have 

MO) = ^ + ^ ^ (T + 1)! r^^-'\ (19) 

So the population of stuck particles decays by an increasingly slow power law as the glass 
transition is approached. Of course, this decay is fastest at early times, so this is when 
the maximum in the mobile peak is formed. The width of the mobile peak grows with the 
usual square-root relation to its mean so, at late times, the population is concentrated in 
two distinct and narrow regions. 

This confirms the result of the above simulation, that the diffusion process in the GR 
model leads, near the glass transition, to a strongly bimodal distribution of displacements. 
Having hopped once, a particle is more likely to hop again in a given time, since relaxations 
in the mobile peak are drawn from the (effectively higher-temperature — Eq. [T^) hop dis- 
tribution \I'hop(T)- Nevertheless, the distribution of residence times throughout the whole 
system (both peaks) is \E'occ(t). Furthermore, the simulational data confirm that, although 
the mean of the mobile peak is a non-trivial function of time and temperature (cf. Eq. |16|) , 
the two peaks together conspire to produce 

N = = t 

(^)hop T 



in agreement with Eqs. |TJ and |T5[ In Ref. p6[, such linearity in time was shown to hold 



in the steady state of any system of random walkers. Nevertheless, in our highly bimodal 
distribution, it appears somewhat remarkable. 
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That the "diffusion coefficient" D should actually comprise an average over two such 
distinct populations might call into question whether the dynamics of the model can properly 
be called diffusion at all: the mean square displacement may be linear in time, but is this 
diffusion in the ordinary sense? To address the question, we calculate the variance in r^(t), 
that is (r^) - (r^)^, and find, for 1 < T < 2, 

- {r'f , /2 , 5r(r) 



(r2) Vrf (3-T)(2-T)v^t^-i' 

The leading term is the standard value for simple diffusion. The way in which this limit is 
approached is temperature dependent, but such sub-dominant terms are unlikely to impinge 
on experimental measurements of the diffusion constant in the ergodic (T > 1) system. The 
deviation from a Gaussian distribution of r at finite times is often characterized by the 



'non-Gaussian parameter' A{t) defined |31| as 

5 (r2) 

in 3 dimensions. Above the glass transition, this is expected [[7| to vanish in the infinite-time 
limit. For 1 < T < 2, we find 

Ai,^ 



(3-r)(2-T) t^-i 

whereas A ~ 1/t above T = 2. So the non-Gaussian parameter in the GR (and Bouchaud's) 
model does indeed vanish, but does so increasingly slowly as the glass transition is ap- 
proached. The same qualitative results were observed by (amongst others) Miyagawa et 
al. in MD simulations ||32[ . 



C. Timescales 

We conclude this section on unforced diffusion in the GR model by summarizing the 
important time-scales that are present in both diffusion and rheology. The diffusion is 
controlled by (T)jjp, equivalent to the mean rate, {t^^)occ^ which vanishes linearly at T = 
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Tg = 1. The mean time between yielding events in the system, (t)^^^ controls the GR 
model's viscosity, and diverges at T = 2, although the mean hopping rate is finite here. The 
median of ^^occIt) also remains finite below T = 2. From Eq. it is the time at which half 
the population has relaxed (when p(0) = |). From Eq. |19| this is 

ti{T) ^[2{T+1)\]T^ (20) 

2 

which has a rather strong divergence at the glass transition (T = 1). Meanwhile, the mode 
of the distributions \l'hop(T) and \1/occ(t), i-e. the most likely time between hops, is always 
zero. 

Given these subtle statistical properties of the self-diffusion process, which arise from 
the power law distribution of relaxation times in the GR model, it is no surprise that the 
GSER (which from Eq. |15| contains only {t)^^^) fails for microscopic probes. 

VI. DIFFUSION IN THE PRESENCE OF SHEAR 

So far we have explored only the linear response functions of the GR model, G*{uj) 
and (r^(t)), and their interrelation. These describe any system for which the applied shear 
stress is either zero or so small as not to disturb thermal equilibrium. Let us now turn 
our attention to applied shear rates which are sufficiently large to infiuence self diffusion, 
and calculate the effect. Of course, particles in Brownian motion are also convected by the 
system's affine shear motion, but this effect (in which we include Taylor dispersion) can 
be subtracted from the overall trajectory to reveal the perturbed diffusive contributionj^. 
(In simple shear, no such subtraction is needed for the component of Brownian motion 



^■^This applies, at least, in numerical studies, and was done in a recent MD simulation |23|. Taylor 
dispersion is the enhanced spreading in the flow direction arising from particles diffusing from one 
streamline to the next. When this and the trivial affine motion is subtracted, the remaining 
Brownian motion had negligible anisotropy. 
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transverse to the flow.) The rheology of the GR model ||TT] , [T2| exhibits strong nonhnearities 
arising from shear-induced hops. These occur as a particle is strained within its well, and 
moves up the energy curve toward the yield value, making a thermal hop far more likely. 
Clearly such shear-induced hops could have a radical effect on the diffusive properties too. 
We assume that, although biased by the affine motion, each such hop still entails a random 
displacement whose statistics is the same as when flow is absent. 

After subtraction of the affine motion, the diffusion constant is defined by 
(^^('^)) t^oo^ 2Dtd in d dimensions. For a random walk of hops, (r^) = Nd so, to 
calculate D, we just require the number of hops executed by an average particle in a given 
long time. That is, 

2D = lim (-) = (21) 



since the time spent in each trap along the walkP^ is drawn from ^E'hopl''"), and each random 
walk becomes typical after a sufficiently long time. We now evaluate Eq. |21| in the presence 
of shear. 

Given a time-dependent macroscopic shear 7(t), the yield rate for a trap which was 
entered at time to was prescribed in section |T| as 

roexp{-/3(E-iA;[7(t)-7(to)]')} 

since the trap has subsequently been sheared by an amount / = •yit) — 7(^0)- At time t, 
the probability of not having hopped from the trap of depth E entered at to (which we call 
fE(t,to)) is therefore 



Mt, to) = exp j-e-^^ / dt' exp (|/3[7(t') - 7(^0)]') 



^^The hop- weighted average {t)^^^, which arises from the microscopic treatment, is completely 
equivalent to the time-weighted average {t^^)^^^^. which arises from the Fokker-Planck equation 
used in Refs. See Eq. |. 
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where Fq = 1 and k = 1 as before. Integrating over the prior distribution of trap depths 
p{E) gives the overall survival probability, 

f{t,to)= (22) 
j^dEp{E) exp [-e-^'^lf exp {^Pb{t') - 7(to)]')} • 

The conditional probability density ipit \ to) of hopping next at time t, given that the present 



trap was entered at to is then 



Ht\to) = -f^f{t,to). (23) 



Hence, the mean yield time, given that the present trap was entered at time to is 

/•OO 

(r),^ = / dt{t-to)^{t\to) 

J to 

f{t,to)dt. (24) 



to 

We next apply Eqs. ^ ^ and to two particular forms of 7(t): steady shear, 7(t) = jt, 
and periodic shear, 7(t) = 7(t + 2n/u). 

A. Steady Shear 

With a time-independent shear rate, the system reaches a steady state. (This is true for 
all temperatures, including those below Tg-. the presence of a finite shear rate destroys weak 



ergodicity breaking |[TT|JT^ .) Therefore, the mean relaxation time becomes independent of 



the time of entry into a trap. Thus Eqs. ^ and ^ are related by 

for any to we might choose. Thus we obtain 

= 2k~^T r du p{-T In ku) ^ (25a) 

JO u 

with J(u) = / rfee""^(^) (25b) 
Jo 

I{0) = f dve"^ (25c) 
Jo 

and = 7V2T (25d) 
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from the substitutions u ^ = kc^^ and 9 = Kt. Here, k is a re-scaled shear rate. 

To proceed further we must approximate J{u) (Eq. |25b| ). Let us define 6q by uI{6q) = 1. 



Then for 6 > 6o, ul{9) grows rapidly with 6*, while for < < it is small. So we can 
write 



Jo Jo 



This approximation is equivalent to saying that shear does not greatly modify the thermal 
activation rate while a particle is in the bottom of a quadratic well, but yielding is immediate 
once the particle has been sheared to the brink of its trap. The result is 

J{u)^r\l/u) (26) 

where is the inverse of I in Eq. |254 We now approximate this inverse function by 



, V for < V < Vo 

/-^(y) - CpU(?/) - r " (27) 

1 + ^/\ny for yo < y 

which is exact in the limits ?/ ^ 1 and y ^ 1. The constant y^ ~ 1.747 solves yo = l + v^ln?/o 
and is the (second]^) point of intersection of the two approximate functions. In fact, it can 
be shown that /approx(l/^) is always greater than J{u) in Eq. |25a| which it approximates. 



So, the resulting approximation for D is a strict lower bound. General expressions for this 
are given in appendix |^. Two regimes of the scaled shear rate arise: n < yo and k > yo- 
Using the exponential prior p{E) we find, in the limit of high shear rate, 

D > — (28) 

So, for fast shear, the diffusion rate is almost proportional to the shear rate. The imposed 
shear drags fluid elements rapidly through trap configurations, with very little time for 



^^Choosing the other solution, yo = 1, creates an upturn in the gradient of I ^, generating 
artefacts. 
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thermal activation. More interesting behaviour occurs in the low-shear regime. With the 
exponential prior, the approximate (lower bound) D is (for T < 3 — see appendix 0) 

D ^ ^^^^ [l - «:^"'c(T)] for «; < y^. (29) 

The function c(T) is given in Fig. |^. 

Below the glass transition (T < 1), the diffusion constant is zero in the absence of shear, 
and particle motion is sub- diffusive, with (r^) growing as some non-integer power of time 
10| . However, any finite steady shear rate ensures a steady rate of yielding events and a 



non-zero diffusion constant. In this regime, from Eq. ^ 

D ^ K^~^. (30) 

"--^ 2Tc{T) ^ ' 

Thus, D is a finite, smooth function throughout all of the (T,7) plane, except for the line 
segment (7 = 0,0<T<1), where it vanishes. For non-zero shear rate, there is no 
singularity at T = 1, where D becomes —1/(2 In k) for low scaled shear rate 



B. Periodic Shear 

A planar periodic shear strain, 7(t) = 7(t+27r/c<j), is experimentally more easily realizable 
than a constant shear rate. If we define the phase to be ut modulo 27t, then we can express 
the mean time (r)^^ spent in each trap on the walk (defined in Eq. ^4]) as a function, r(0), 
only of the phase at which the current trap was entered. This is done in appendix ^ by 
averaging (r)^^ over all oscillatory cycles, for all times to corresponding to phase (p. If p(0) 
is the probability that a trap was entered at phase of the cycle, then the mean residence 
time in Eq. ^ is 

(^)hop= r#p(</')r(0). (31) 

Let p{(j) I 0o) be the conditional probability density of a hop occurring at phase 0, given that 
the previous hop was at 0o- Then p{(j)) is the (normalized) solution to the integral equation 
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p(0)= / P(0|0o)p(0o)#o. (32) 

JO 

The quantities p(0 | 0o) and t(0) have simple relations to the distribution ilj{t \ to) given in 
Eq. In appendix ^ they are given in terms of a single cycle of the imposed waveform. 

Equation ^ is a homogeneous Fredholm equation of the second kind, for which the 
general solution is not known. However, the problem greatly simplifies in the limit of small 
amplitude oscillations, which we now treat. (Note that this perturbative treatment is only 
possible for T > 1, where an unperturbed steady state exists.) In the absence of shear p(0) 
and t(0) must be independent of 0. In that case, we write p{(f)) = Xjln and t(0) = tq. If 
the oscillatory shear creates only small perturbations Ap(0) = p(0) — l/27r and Ar(0) = 
r(0) — To, then Eq. |3l|, to first order in small quantities, becomes 



^rr{^)d<f>. (33) 
Zrr Jo 



We have used the fact that p is normalized, so that J^^ ^p{4>) d(j) = 0. Hence we need only 
expand Eq. ^to lowest order in strain amphtude and average it over a cycle to find {t)^^^ 
and hence the diffusion constant. The result for an exponential prior p{E) = exp{—E) and 
sinusoidal shear of amplitude 70 is 

^-..JI__,|.-^"-*^*^OhS (34) 

For T < 3, the integral is convergent as — * 0. In fact, for any waveform, the low-frequency 
limit for an exponential prior is 

^ ^ - I'oU^'^-'giT) + O(7o^). (35) 

where, in the sinusoidal case, g(T) = (vr/2) cosec [(T — l)7r/2]. In the high-frequency limit, 
Eq. ^ yields instead 

For larger amplitudes, we can use the results of section |V1 A| . Given a sufficiently low 
frequency {uj D), the shear rate at any given point on the cycle is almost constant. In 
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that case, the steady-shear expressions for D ^ can be averaged over a sinusoidal cycle, 
yielding (for u <ti D) 

' ■ ^l^'' < 2T (36) 

and, to within logarithmic factors, 

Docujo ; 7^cj2>2r (37) 

where h{T) = 2^^-^^/^c{T) 5(T/2, T/2)/7r, c is given in Fig. | and B is Euler's integral of 
the second kind. 

We now summarize the effects, calculated above, of low-frequency periodic shear on 
diffusion just above the glass transition. For the smallest amplitudes 70, we expect the 
change in the diffusion constant to be proportional to 'JqUj'^^^. At larger amplitude, this 
should cross over to (700;)"^"^, and at the largest amplitudes, D is almost proportional 
to uj'jo- We interpret these three regimes as follows. The smallest oscillations serve only 
to vibrate particles within the bottoms of their quadratic potential wells, thus effectively 
making all wells just a little shallower. Larger amplitude strains cause yielding of even the 
deepest wells in the time taken for the affine shear to reach the yield point. However, the 
part of the population in the shallowest wells are thermally activated more quickly than this 
affine shear time-scale. The dividing line between those two parts of the population depends 
on the shear rate. At the highest shear rates, the global affine shear accounts for almost all 
yielding. 
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VII. DISCUSSION AND CONCLUSION 

We have studied quantitatively, within a simple hopping model of glassy dynamics, the 
relation between diffusion and rheological responses (in the linear response regime) and (for 
the nonlinear regime) the coupling between these two aspects of the dynamics. 

In the presence of a broad distribution of relaxation times, as the model possesses (and 
as is generic in hopping models of the glass transition) one can expect strong violations of 
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the generalized Stokes-Einstein relation (Eq. P when applied [0] to the self-diffusion of 



representative particles in the medium (small probes). The relation holds only for probe 
particles large enough that their surroundings are properly viewed as a continuum [p!^ . 



Effective-medium theories, and also simple forms of mode-coupling theory (in which 
a single mode of slow relaxation dominates) can give misleading predictions under these 
conditions. For glassy systems, use of the GSER will under-predict the diffusivity of particles 
and/or the rheological viscoelastic moduli. This is because the latter are dominated by the 
most immobile and the former is dominated by the most mobile particles. 

In fact, the nature of diffusion in the model is quite subtle (see section close to 
(but above) the glass transition one has an apparently bimodal behaviour. Starting from 
any initial equilibrium state, a fraction of particles remain stuck for a long time, but any 
particle that has hopped once remains mobile thereafter. It is peculiar that, despite this, the 
diffusive behaviour conspires to be relatively normal (at least for low order moments of the 
displacement distribution; moments of fractional negative order would presumably reveal a 
different story). 

An additional peculiarity of diffusion, which our model does not include, but which could 
also lead to enhanced diffusivity near the glass temperature, has been observed in a recent 
molecular dynamics simulation |3^. In Ref. the effective dimensionality of the most 



mobile random walks was observed to decrease with temperature. At low temperature, 
particles moved along string-like clusters (also seen in Ref. [^), and therefore covered 
greater distances than in an uncorrelated random walk in three dimensions. 

We have also applied our GR model to calculate the change in the self- diffusion constant 
when a material is sheared. After subtraction of affine motion (including Taylor dispersion) 
23| , or equivalently restricting attention to diffusion perpendicular to the shear direction. 



one finds a strong effect of imposed flow on the mean jump rate [|I^ and hence on the diffusion 
constant. The effect is particularly extreme below the glass transition, where an anomalous 
(sub- diffusive) behaviour is converted to a finite diffusivity which has, instead, a power 
law dependence on the steady shear rate. Such effects should be accessible in scattering 
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experiments (with wavevector almost perpendicular to the flow direction) on labelled small 
probes. In oscillatory shear, an enhancement in D was predicted in section |V| for systems 
above (but near) the glass temperature; this could also be probed via scattering. 

The GR model as described in section |IIT| was originally developed in Ref. [|11|] to repro- 
duce the generic rheology of a class of "soft glassy materials" . (We comment further below on 
its relevance to conventional glasses.) This class was argued to include, for example, foams, 
emulsions, pastes and slurries. Experimentally, their linear viscoelastic behaviour is often 
characterized by a nearly constant ratio of the elastic and loss moduli G"(c<j), G"{uj) {G" /G' 
is usually about 0.1) with a frequency dependence that is either a weak power law (clay 
slurries, paints, microgels) or negligible (tomato paste, dense emulsions, dense multi-layer 
vesicles) . This behaviour persists down to the lowest experimentally accessible frequen- 
cies. Sometimes a regime is seen at small u where G' is constant and G" is decreasing (which 



can be interpreted in terms of the model's behaviour below Tg |11,12,2 



As mentioned in sections | and [TTII , there are two differences between the GR model used 



in this paper and its soft counterpart in |TT[]. The first is that the soft GR model refers not to 
particles, as we have done, but to mesoscopic material "elements" , large enough for a local 
elastic strain variable to be defined but small enough to have strong heterogeneity in local 
yield energies. (For the case of a foam, say, an "element" could correspond to a domain of 
several bubbles, and "yield" to a local topological change.) It is not necessarily clear what is 
meant by self-diffusion of such elements, so in the present paper we have retained a particle 
picture, although this is more natural for conventional glass-forming liquids than for soft 
glassy materials. 

The second difference is that in the soft GR model, T is replaced by an effective "noise 



temperature" x. It was argued in |TT|] that the resemblance to thermal activation is formal: 
the "activated" yield processes are viewed as arising primarily by coupling to structural 
rearrangements elsewhere in the system. Indeed, the elastic energies associated with local 
rearrangements in foams and the like are many orders of magnitude in excess of ksT, so 
any interpretation of x as a true temperature is somewhat unconvincing for these materials 
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(in contrast to conventional glasses). Although this interpretation remains problematical, as 



discussed in [|ri|,|T^ , with it the soft GR model is able to reproduce many of the rheological 
properties of soft glassy materials. 

The results we have obtained for the GR model, concerning breakdown of the GSER and 
the nature of shear-induced diffusion, equally apply (with T ^ x) to the soft GR model. 
The breakdown of GSER, arising from the fact that diffusion and rheology probe different 
aspects of the relaxation spectrum, is equally natural in this case; the other contribution 
to its breakdown discussed in [IV A| , arising from the temperature dependence of the local 



attempt frequency, even more so (since x is anyway not a true temperature). The results of 
section ^ for shear-induced diffusion could also be quite interesting for soft glassy materials, 
in which it is easy to apply shear strains large enough to strongly perturb the intrinsic 
relaxation times. The length scales in these materials can be probed via light scattering; 
index matching is often possible so that true tracer diffusion (of a small subset of unmatched 
droplets or particles) can be measured For oscillatory flows, an important innovation 
is the echo technique in which the positions of a given scatter at identical points in the 
shear cycle are compared. 

Results on dense emulsions PB| suggest, in fact, that not only are shear-induced reor- 



ganizations easily detectable, but that these have strong temporal and spatial correlations 
— regions that reorganize at a given point in the shear cycle will do so again in the next 
one. Such correlations are not included in the (soft) GR model(s) although they might be 



added in principle |Q. (It would require escape from a shallow trap to be preferentially 
into another shallow trap, though we do not advocate appending such ad-hoc correlations 
to this simplified model.) It would be very interesting to know whether the same applies 
in conventional glasses; preliminary work on colloidal suspensions (which are traditionally 
thought of in these terms) suggests not [ p^ . 

Finally, we return to the phenomenology of the GR model in relation to conventional 
glass-forming liquids. As explained in Section [Lll[B, the GR model in its basic form predicts a 



viscosity divergence at T = 2, while the diffusion constant vanishes only at T = 1, where the 

31 



system has a glass transition to a non-ergodic state. This existence of a temperature range 
with finite diffusivity but infinite viscosity appears to be at odds with the experimentally 
observed behaviour of conventional glasses. A common divergence of viscosity and inverse 
diffusion constant can however be incorporated into the GR model through a cutoff -Emax 
on the distribution of yield energies p{E); this modification of the model has in fact already 
been discussed in Ref . . It yields a viscosity which initially follows the original power-law 
divergence as T = 2 is approached, but then crosses over to ~ exp{E^s.x/T) as T is lowered 
further. Similarly, the predicted 1/D would first seem to diverge as T = 1 is approached 
from above, but actually remain finite there and eventually approach infinity at the same 
temperature as the viscosity (T = 0). The shear moduli obey the power laws (Eq. ^ down 
to a cutoff frequency cUmin = ^w{~Ema.x/T), but then cross over to low-frequency Maxwell 
behaviour (G" ~ cu^, G" ~ u). 

The above simple temperature dependences apply if, as we did throughout, we assume 
that the (prior) density of yield energies is temperature independent. This is of course an 
approximation; Odagaki, for example, suggested that the width of p{E) may in fact scale as 
the inverse of the amount of free volume v{T) in the system In all temperature depen- 
dences, Tv(T) then replaces T. As a consequence, if the free volume decreases smoothly 
to zero at a finite temperature Tvf, a Vogel-Fulcher-like divergence of rj and 1/D at Tvf 
(rather than the above Arrhenius behaviour) would be predicted by the GR model with 
energy cutoff. One would then be inclined to locate the glass transition at that point; if 
we revert to temperature independent p{E), this corresponds to T = (rather than T = 1, 
which is the appropriate choice in the absence of an energy cutoff). Results in the range 
< T < 1 for the GR model with cutoff may therefore actually apply to supercooled liquids 
above the glass transition. These include a dynamic modulus G*{u!) which becomes more 
Maxwellian in the low- frequency range as T is lowered (for T < 1, one has G' ~ const, and 
G" ~ uj'^^^ above the cutoff frequency tUmin [0, and hence Maxwell behaviour with relax- 
ation time l/cUmin for T — 0), in an intriguing correspondence with data taken by Menon 
et al. 
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APPENDIX A: GENERAL EXPRESSIONS FOR THE DIFFUSION CONSTANT 

UNDER STEADY SHEAR 



Substituting Eqs. ^ and into p5a| yields lower bounds on D (upper bounds on D ^) 
in two regimes: 



--1 







< 2n-'T I du ^"""^ (l + Vh[^] 

for K > Ho (Al) 



Pi- 


-T In nu) 
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Pi- 


-T In ku) 
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Pi- 


-T In nu) 


V? 



for K < uq. (A2) 



In each case, D is approximately equal to the lower bound for T < 3. 

APPENDIX B: EVALUATION OF THE PHASE-DEPENDENT RELAXATION 
TIME AND CONDITIONAL PROBABILITY UNDER PERIODIC STRAIN 

The probability of leaving a trap at time t given that it was entered at to (at phase 
00 = i^^o mod 27r of the oscillatory cycle) is '0(t | to)- Hence the mean time spent in the trap, 
given that it was entered at phase 0o; is 



which, from Eq. p^, gives 



^('/'o) = / it-<Po/^) ^(i|0o/t^) dt 



r(0o) = / f{t,<Po/uj)dt 



with / defined in Eq. For compact notation, we define a{ujt) = 7(t)/V2T and te(0o) 
according to 
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r(0o) = / p{E)TE{(l)o)dE. 
Jo 



Thus we find 



(Bl) 



dc/) exp <— M / d(j)' exp[a(0') — a(0o)]^ > 

'0 [ J 



Since exp[a(</)') — a(0o)]^ is a periodic function of (p, the integrand in Eq. ^T] increases by 
a constant factor with each cycle — > + 2tt. Thus the R.H.S. is a geometric series of 
integrals over a single cycle, 



OO r / »27r 

V exp -u / ci0' exp[a(0') - a{(f)o)]' 

Z^al \ Jo 



n=0 

(t>0 + 2TT 
<t>0 



exp < —u i d(j)' exp[a(0') — a(0o)]' 



Finally, r(0o) is expressed in terms of a single cycle as 

T /•^-' p(-T In cum) /•2; 



LJ Jo 



du 



u 



X 



(B2) 



1 — exp |— M ^Q^dcj)' exp[a(0') — a(0o)]^} 

where 6(0o ~ 0) is the Heaviside step function. 
Similarly, 



p(0 I 0o) = / 5([co'tmod27r] — 0) '0(t | 0o/ti;) (it 

(with 5{x) the Dirac delta function) is also a geometric series, from which it follows 

d 



p(0 I 0o) = ^ / dup{—T\nuju) 
Jo 



X 



(B3) 



exp 



(/>+27r0((/-o-0), 

4>o 



'rf0' exp[a(0') -a(0o)]^} 



1 - exp {-u S^^dct>' exp[a(0O - a{<t>o)Y} 
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FIG. 1. Normalized distributions of number of hops performed, pt{N), in simulations at various 
temperatures T and times t. For (T = 1.1, t = 1000), pt{0) = 0.476 ± 0.002 and the rest of the 
distribution is magnified in the inset. 
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FIG. 2. The function appearing in Eqs. 29 and^, c(T) = T^^^q^^— ^T^2 (T— l)erfc-v/TTnyO) 
where erfc is the complementary error function and yo ^ 1-747. 
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